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(§) Methode de simulation pour predire en fonction du 
temps une composition detailtee a Q composants et/ ou 
pseudo-composants d'un fluide produit par un reservoir. 
Elle est du type consistant a: 

- representer le reservoir sous la forme cf un reseau de 
mailles (j) dont chacune constitue un volume elementaire 
rempli de fluide, 

- realiser, de facon conn ue en so i, une simulation com- 
positionnelle du fluide regroupe a N composants et pseudo- 
composants (i), N etant inferieur a Q, et elle est caracterisee 
en ce qu'eile consiste en outre a: 

e) determiner, a chaque pas le temps (m) et pourcha- 
que maille (j) la composition des phases liquide et vapeur du 
fluide detaille de la composition Q du fluide, a partir des va- 
leurs de la fraction vaporisee (6j m ) et des constantes d'equi- 
Itbre (k t f) du fluide regroupe, 

f) evaluer pour chaque maille, au pas de temps (m+1), la 
quantite molaire de chacun des Q composants et/ ou pseu- 
do-composants du fluide detaille a partir des valeurs corres- 
pondantes au pas de temps (m), des debits des phases du 
fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a I'etape e), et a 

g) evaluer, pour chaque puits de production, la compo- 



sition detaillee du fluide produit entre des instants t et t* cor- 
respondant aux pas de temps m et m+1 , a partir des debits 
des phases du fluide regroupe et des compositions des pha- 
ses liquide et vapeur du fluide detaille determinees a I'etape 
e). 
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Methode de simulation pour predire en fonction du temps 
une composition detaillee d'un fluide produit par un reservoir 

5 

La presente invention concerne une methode de simulation pour 
predire, en fonction du temps, la composition detaillee dun fluide produit 
par un reservoir et plus particulierement la composition detaillee d'un fluide 
contenu dans et produit par un gisement petrolier dans lequel sont implantes 

10 un ou plusieurs puits de production. 

La simulation compositionnelle d'un gisement petrolier est 
couramment utilisee pour fournir les profils previsionnels de production du 
gisement qui permettent notamment de determiner le schema de production 
le mieux adapte a ce gisement. 

15 La simulation compositionnelle d'un gisement petrolier est mise 

en oeuvre non pas en utilisant une description reelle du fluide mais en 
utilisant un fluide modelise par un nombre de composants plus reduit que le 
nombre de composants du fluide reel. En effet, le nombre de composants du 
fluide de gisement etant relativement grand, la modelisation avec tous les 

2 0 composants conduirait a des temps de calcul trop importants. Ce temps de 

calcul prohibitif a amene les specialistes a regrouper les composants purs en 
pseudo-composants, par exemple des pseudo-composants regroupant Tazote 
et le methane, des pseudo-composants regroupant les hydrocarbures en C3 
et C 4 , ... etc., et a realiser la simulation de gisement sur un nombre reduit 
25 N de composants purs et de pseudo-composants. Une telle composition 
reduite est designee par Texpression "composition regroupee du fluide" ou 
"fluide regroupe". Generalement N represente 5 a 8 composants et pseudo- 
composants, ce qui est considere comme etant suffisant pour bien 
representer le comportement du fluide de gisement aux conditions de fond. 

3 0 Avant d'effectuer la simulation compositionnelle de gisement, on 

represente le reservoir ou gisement sous la forme d'un reseau de mailles 
dont chacune constitue un volume elementaire dudit reservoir. Le nombre 
de mailles peut atteindre plusieurs milliers et chaque maille presente des 
proprietes qtii lui sont propres, comme la geometrie, la porosite, la 
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permeabilite, etc... De plus, a certaines mailles peut correspondre au moins 
un puits d'injection ou de production qui est implante dans ledit reservoir. 

La simulation de gisement permet de calculer pour chaque maille 
un certain nombre de variables principales du fluide de gisement, lesdites 
variables pouvant etre la quantite (nombre de moles) de fluide, la compo- 
sition ou fraction molaire de chaque composant et pseudo-composant du 
fluide regroupe et la pression regnant dans chaque maille. Ces variables 
sont connues a l'instant t=0 (debut de Texploitation du gisement) et sont 
ensuite calculees par la simulation de gisement a chaque instant t ou pas de 
temps m. Pour chaque maille, il est possible, a partir de ces variables 
principales, de calculer toute autre propriete du fluide present dans la 
maille, comme le nombre de phases, la composition de chaque phase, etc. 
Pour chaque maille a laquelle est associe un puits de production, il est 
possible, a partir de ces variables principales et des contraintes de 
production imposees, de calculer egalement le debit de production et la 
composition du fluide produit par ledit puits. 

Les proprietes thermodynamiques du fluide peuvent etre 
calculees en utilisant une des equations d'etat bien connues des speciaJistes 
et qui ne seront done pas rappelees. 

La modelisation a N composants, si elle est suffisante pour 
representer le comportement du fluide de gisement aux conditions de fond, 
n'est plus appropriee pour simuler le comportement du fluide dans les 
installations pour son exploitation en surface, laquelle exploitation necessite 
de connaitre une composition plus detaillee a Q composants et/ou pseudo- 
composants du fluide produit par le gisement, Q etant superieur a N et par 
exemple de Tordre de 20 a 30. 

Jusqu'ici, on obtenait cette composition detaillee a Q composants 
et/ou pseudo-composants a partir de la composition regroupee N en 
supposant que la composition des pseudo-composants restait constante dans 
le temps. 

Une telle maniere de proceder est source d'erreurs car la 
composition des pseudo-composants varie dans le temps, notamment en 
fonction de la pression du gisement ou suite a l'injection d'tm gaz dans le 
gisement. 
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De ce fait, la composition detaillee obtenue reste tres approxima- 
tive et ne permet pas de predire de maniere satisfaisante le comportement du 
fluide dans les conditions Sexploitation en surface bien des annees apres le 
debut de la production. 
5 La presente invention a pour but de proposer une methode de 

simulation pour predire la composition detaillee du fluide produit par un 
reservoir en tenant compte des parametres du reservoir et qui soit beaucoup 
plus precise. 

^invention a pour objet une methode du type consistant a : 
10 a) representer le reservoir sous la forme d*un reseau de mailles (j) 

dont chacune constitue un volume elementaire rempli de fluide, 

b) definir le fluide par une modelisation regroupee a N 
composants et pseudo-composants (i) et a determiner une equation d'etat 
decrivant le fluide dans cette modelisation regroupee, 
15 c) definir egalement le fluide par une modelisation detaillee a Q 

composants et/ou pseudo-composants, Q etant superieur a N, et a deter- 
miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa9on connue en soi, une simulation composi- 
2 0 tionnelle du fluide regroupe a N composants et pseudo-composants (i), 

ladite simulation compositionnelle permettant de calculer au moins pour 
chaque maille (j) et a des pas de temps consecutifs (m, m+1, ...) la fraction 
vaporisee (©™), les constantes d'equilibre liquide-vapeur (k™) de chaque 

composant (i), les debits d'injection ou de production (s™) et pour chaque 
25 paire de mailles (j, h) les debits des phases liquide (u£ h ) et vapeur 

fluide a N composants et pseudo-composants, et elle est caracterisee en ce 
qu f elle consiste en outre a: 

e) determiner, a chaque pas de temps (m) et pour chaque maille 
(j) la composition des phases liquide et vapeur du fluide detaille de la 

30 composition Q du fluide, a partir des valeurs de la fraction vaporisee (0™) et 

des constantes d'equilibre (k™) du fluide regroupe, 

f) evaluer pour chaque maille, au pas de temps (m+1), la quantite 
molaire de chacun des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (m), des debits 
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des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a Tetape e), et a 

g) evaluer, pour chaque puits de production, la composition 
detaillee du fluide produit entre des instants t et t 1 correspondant aux pas de 
temps m et m+1, a partir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
1'etape e). 

Selon une autre caracteristique de la presente invention, les 
etapes e, f et g sont mises en oeuvre en meme temps que la realisation de la 
simulation compositionnelle du fluide regroupe. 

Selon une autre caracteristique de la presente invention, les 
resultats de 1'etape d sont memorises puis utilises ulterieurement pour la 
mise en oeuvre des etapes e, f et g. 

La methode selon la presente invention est mise en oeuvre pour 
un reservoir constitue, par exemple, par un gisement petrolier. De maniere 
connue, on realise une simulation de gisement. A cet effet, on represente le 
gisement sous la forme d'un reseau maille, certaines des mailles ou un 
groupe de mailles etant associes a un des puits de production implantes dans 
le gisement petrolier a exploiter. 

Pour reduire le temps de calcul, la simulation compositionnelle 
du gisement est effectuee sur un nombre limite N de composants et pseudo- 
composants, par exemple de 5 a 8, definis de la maniere indiquee 
precedemment, ces composants et pseudo-composants etant selectionnes en 
fonction de la nature du gisement. 

Le fluide regroupe a N composants et pseudo-composants est 
decrit par une equation d'etat qui pourrait etre par exemple celle de PENG- 
ROBrNSON adaptee a ladite composition regroupee. 

Au debut de I'exploitation du gisement, on connait, par des 
mesures effectuees au prealable, pour chaque maille j du reseau, les 
variables principals a Tmstant t=0 ou pas de temps m=0. Ces variables 
principales sont : 
la quantite de fluide : f? 

la fraction molaire du composant i : avec 1 < i < N 
la pression : pj" 
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Pour passer du pas de temps m au pas de temps m+1, la 
simulation compositionnelle comprend plusieurs etapes detaillees ci-apres. 

Dans une premiere etape et pour le pas de temps m, un calcul 
d'equilibre liquide-vapeur (flash) est effectue sur le fluide regroupe, pour 
5 chaque maille j afin de determiner : 

la fraction vaporisee 0™ 

la fraction molaire de chaque composant i dans la phase huile (x™ avec i = 1...n) et 
la fraction molaire dudit composant i dans la phase gaz (y™ avec i = 1...n). 

Ces diverses fractions sont determinees avec le systeme 
1 0 d'equations suivant : 

O^x-^.-.x-Jet ^(y^-.y^avec l<i<N 

et 

N N 

1=1 i=i 

15 

dans lesquelles : 

<3>j est la fugacite du composant i. 

Selon une premiere caracteristique de l'invention, on memorise la 

y m 

fraction vaporisee 0™ et les constantes d'equilibre k™ = — ^ avec i variant 

20 entre 1 et N, afin de les utiliser dans une autre etape denommee 
"delumping" qui sera explicitee ulterieurement 

Dans une deuxieme etape et pour chaque maille, on evalue 
certaines des proprietes de chaque phase liquide et gaz, a savoir masse 
volumique (p™ et p™), viscosite, saturation et permeabilite relative. 

25 Dans une troisieme etape, pour chaque maille j et au pas de 

temps m, on calcule les coefficients de l'equation pour la pression. 

Les equations pour la pression d'une maille j a 1'instant t+At ou au 
pas de temps m+1, sont : 
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20 



c 4 P '\+ P ' r Zn(pr'-pr')+sr 

l At } h€j(j) 



Pour resoudre ces equations, les coefficients a calculer sont : 
5 - la compressibilite totale (fluide + roche), 

- les transmissivites generalisees 

TS = £a?(x o x j+ Jl B yJ m tjh 

N 

10 - le teime source ; s? = V s ro 

Dans ces equations : 
h sont les mailles voisines de la maille j, h e J(j), 

a? est le volume molaire partiel du composant i dans la maille j au pas de 
temps m, X™ et \™ sont les mobilites des phases huile (o) et gaz (g), 
15 t jh est la transmissivite entre les centres des mailles j et h, 1 etant la maille 
amont. 

Dans une quatrieme etape, on resout pour toutes les mailles les 
equations pour la pression a Tinstant t+At c'est-a-dire au pas de temps m+L 
Ces equations forment un systeme lineaire du type ; 



A. p = B 



oup = (p, m * l ,...p™* , ) T , J etant le nombre de mailles du reseau. 

Ces equations sont resolues par un procede iteratif standard. 
25 Dans une cinquieme etape, on calcule au pas de temps m+1, les 

debits des phases entre chaque paire de mailles adjacentes et pour chaque 
puits associe a partir, par exemple, de la loi de DARCY. 

U ojb --Poj A.ol -t jh ^p h ~p i ) 

U oh --Pa\ tihlPh "Pj ) 

30 
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1 etant la maille amont avec 1-j ou h. 

Selon une autre caracteristique de l'invention, on stocke ou 
memorise les debits des phases qui seront utilises dans le delumping. 

Dans une sixieme etape, on evalue la quantite f et la composition 
z du fluide regroupe, dans chaque maille j et au pas de temps m+1 a partir 
des valeurs au pas de temps m, des debits des phases evalues a la cinquieme 
etape et des compositions des phases calculees a Fetape 1 . Pour ce faire, on 
applique les lois de conservation definie par les equations suivantes : 



f-' = f» - At] £ (u™ + u~) - sA 

|>:Hj) J 



10 Par un processus iteratif, on repete les etapes ci-dessus pour 

calculer les variables principales aux temps t+2At, t-^3At, etc... ou encore 
aux pas de temps m+1, m+2, m+3, etc., rincrement At pouvant etre 
constant ou variable. 

Selon la presente invention, on definit une modelisation detaillee 
15 a Q composants et/ou pseudo-composants ainsi qu'une autre equation d'etat 
decrivant le fluide ainsi modelise. 

Pour cela, on choisit prealablement les composants et/ou pseudo- 
composants dont on veut connaitre les caracteristiques et qui sont 
necessaires pour representer le fluide de production dans les conditions de 

2 0 surface. Le nombre Q qui est superieur au nombre N de la composition 

regroupee est generalement de Tordre de 20 a 30. On connait la composition 
du fluide detaille a Q composants et/ou pseudo-composants au pas de temps 
m=0. 

Les donnees stockees aux premiere et cinquieme etapes sont 
25 utilisees pour les etendre aux Q composants et/ou pseudo-composants de la 
composition detaillee. 

Conune Tecoulement peut s'effectuer dans n'importe quelle 
direction entre les mailles voisines, chaque maille ne peut pas etre traitee 
independamment. C'est pourquoi, a chaque pas de temps, toutes les mailles 

3 0 sont traitees ensemble. 
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Pour realiser le delumping, on fait lltypothese qu'a chaque pas de 

temps m, 

- la fraction molaire de la phase vaporisee dans chaque maille est 
independante du nombre de constituants du fluide, 



0™ = e™ 



le nombre de moles de liquide et de vapeur s'ecoulant entre les mailles et 
dans ou hors des puits est independant du nombre de constituants du 
fluide, 



U- h = u™ et U^h = u m 



Les lettres en majuscule sont utilisees pour representer le fluide 
detaille et les lettres en minuscule sont utilisees pour representer le fluide 
regroupe. 

Le calcul des constantes d'equilibre peut etre effectue a partir de 
1'expression ; 

Ln[K, ] = AC 0 + ± AC n . n m ou AC n = Cj; - C n v 

comme indique par C. LEIBOVICI dans son article "A Consistent 
Procedure for Pseudo-Component Delumping" paru dans FLUID PHASE 
EQUILIBRIA, 117(1996), 225-232. 

Pour une equation d'etat a deux parametres, comme celle de 
PENG ROBINSON, l'equation generalisee de LEIBOVICI peut etre ecrite 
pour le fluide detaille de la maniere suivante : 

MK I 7) = i 1 ?+p-^+yr b i 

ou ai et bj sont les parametres de l'equation d'etat pour le composant I et ou 
les parametres r£\p™ ety™ sont calcules a partir des constantes d'equilibre 

pour le fluide regroupe (k™) en minimisant la fonction : 
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De ce fait, lorsqu'on connait la composition detaillee Zi7 dans la 
maille j au pas de temps m, on peut estimer les compositions detaillees des 
phases liquide et vapeur dans la meme maille j et pour le meme pas de 
temps m. 

5 

XI? = ^ r Y™ = M r avec 1 < I < Z 

Ij l + ©7(K"-l) ' l + G^K^-l) 

I etant le composant et/ou pseudo-composant de la composition 
detaillee. Ces compositions peuvent etre normalisees afm que leur somme 
10 soit egale a 1 si necessaire. 

Lorsqu'on connait les compositions des phases detaillees pour 
toutes les mailles du reseau au pas de temps m, il est possible d'estimer les 
compositions detaillees zip 1 au pas de temps consecutif m+1 a partir des 

equations ci-apres : 



15 



z$?Fr-(YSS3+x?s£)- 2(Yi?U5. + X!?USh) 
zr 1 = Z^r* avecl= 1,...,Q 



ou 

20 ft' = f?-(s;+sj)- E(ua + u3.) 

hcJ(j) 

dans lesquelles j'=j pour uri debit positif de la maille j vers la maille h ou 
dans le puits et j'=h pour un debit negatif, j* correspondant au fluide injecte 
pour les puits d'injection ou S sera negatif. 

Grace a la presente invention, on obtient une composition 

25 detaillee a Q composants et/ou pseudo-composants qui permet de mieux 
definir les profils de production previsionnels et ce, avec une precision 
accrue comparee aux methodes anterieures. Cette plus grande precision des 
profils de production previsionnels permet un choix de schema de 
developpement (type et taille d'installations de surface, nombre de puits, 

30 etc..) plus fiable, avec pour consequence une augmentation de la rentabilite 
et une diminution du risque economique. 
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Dans ce qui precede, certaines valeurs obtenues au cours de la 
simulation de gisement sont stockees ou memorisees pour etre utilisees 
ulterieurement dans le dehimping. II va de soi que ces memes valeurs 
pourraient etre utilisees directement si on decide d'effectuer le delumping en 
meme temps que ladite simulation de gisement. 
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REVINDICATIONS 



5 1. Methode de simulation pour predire, en fonction du temps, la 

composition detaillee d'un fluide produit par un reservoir, du type consistant 

a : 

a) representer le reservoir sous la forme d'un reseau de mailles (j) 
dont chacune constitue un volume elementaire rempli de fluide, 
10 b) definir le fluide par une modelisation regroupee a N compo- 

sants et pseudo-composants (i) et a determiner une equation d'etat decrivant 
le fluide dans cette modelisation regroupee, 

c) definir egalement le fluide par une modelisation detaillee a Q 
composants et/ou pseudo-composants, Q etant superieur a N, et a deter- 

15 miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa?on connue en soi, une simulation composi- 
tionnelle du fluide regroupe a N composants et pseudo-composants (i), 
ladite simulation compositionnelle permettant de calculer au moins pour 

2 0 chaque maille (j) et a des pas de temps consecutifs (m, m+l, ...) la fraction 
vaporise {q™) 9 les constantes d'equilibre liquide-vapeur (k^ 1 ) de chaque 

composant (i), les debits d'injection ou de production (s™) et pour chaque 

paire de mailles (j, h) les debits des phases liquide (u^ h ) et vapeur (u" h ) du 

fluide a N composants et pseudo-composants, et elle est caracterisee en ce 

2 5 qu'elle consiste en outre a : 

e) deteiminer, a chaque pas le temps (m) et pour chaque maille 
(j) la composition des phases liquide et vapeur du fluide detail le de la 
composition Q du fluide, apartir des valeurs de la fraction vaporisee (dj) et 

des constantes d'equilibre (k™) du fluide regroupe, 

3 0 f) evaluer pour chaque maille, au pas de temps (m+1), la quantite 

molaire de chacun des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (m), des debits 
des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a 1'etape e), et a 
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g) evaluer, pour chaque puits de production, la composition 
detaillee du fluide produit entre des instants t et t' correspondant aux pas de 
temps m et m+1, a partir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
5 l'etape e). 

2. Methode selon la revendication 1, caracterisee en ce que les 
etapes e, f et g sont mises en oeuvre en meme temps que la realisation de la 
simulation compositionnelle du fluide regroupe. 

3. Methode selon la revendication 1, caracterisee en ce que les 
10 resultats de Tetape d sont memorises puis utilises ulterieurement pour la 

mise en oeuvre des etapes e, f et g. 



rfSDOCID: <FR 2775094A1_1_> 



2775094 



t 



REPUBLIQUE FRANQAISE 



1NSTITUT NATIONAL 
de la 

PROPRIETE INDUSTRIELLE 



RAPPORT DE RECHERCHE 

PRELIMINA1RE 

etabli sur la base des demieres revendications 
deposees avant le commencement de la recherche 



N* d*en regis trement 
national 



FA 556194 
FR 9801970 



Categorie 



DOCUMENTS CONSIDERES COMME PERTINENTS 



D,A 



Citation du document avec indication, en cas de besom, 
des parties pertinent es 



C. LEIBOVICI ET AL: "A Consistent 
Procedure for Pseudo-Component Delumping" 
FLUID PHASE EQUILIBRIA. 

vol. 117, 1996, pages 225-232, XP002083980 

* le document en en tier * 

J. A. TRANGENSTEIN ET AL: "Mathematical 
Structure of Compositional Reservoir 
Simulation" 

SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL 
COMPUTING, 

vol. 10, no. 5, septembre 1989, pages 
817-845, XP002083981 

* le document en entier * 



Revendications 
concents 
de la demande 
examnee 



1-3 



OOMAINES TECHNIQUES 
RECHEHCHES (lnl.CL.6) 



G06F 



Date d'achevement de la recherche 

12 novembre 1998 



Examinateur 

Abram, R 



CATEGORIE DES DOCUMENTS CITES 

X : particulierement pertinent a tui seul 

Y : particuJierement pertinent en comWnarson avecun 

autre document de la meme catdgone 
A : pertinent a rencontre d"au moins une revindication 

ou arr t ere- plan tecrtnologicpje general 
O : divulgation non-ecrite 

P : document intercalate 



T : theone ou princtpe a la base de I invention . 

E r document de brevet beneticiant dune date anterieure 

a ta date de depot et qui n'a ete publiequa cette date 

de depot ou qu a une date posterieure. 
O : cite dans la demande 
t : cite pour cfautres raisons 

& : membre de ta meme famine, document oorrespondant 



BNSDOCID: <FR 2775094A1_I_> 



